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ABSTRACT 

We show that the gravitational magnification factor averaged over ah configurations 
of lenses in a locally inhomogeneous universe satisfy a second order differential equation 
with redshift z by taking the continuous limit of multi-plane gravitational lens equation 
(the number A'l of lenses oo) and that the gravitationally magnified Dyer-Roeder 
distance in a clumpy universe becomes to that of the Friedmann-Lemaitre universe for 
arbitrary values of the density parameter f^o and of a mass fraction a (smoothness 
parameter) . 

Subject headings: cosmologyitheory — distance scale — gravitational lensing 



1. Introduction 

A light ray propagation in a locally inhomogeneous universe has been investigated by many 
authors using analytical and/or numerical methods (e.g., Omote & Yoshida 1990; Yoshida & Omote 
1992; Schneider, Ehlers & Falco 1992, and references therein). By taking gravitational lens effects 
into account, Weinberg (1976) showed that in a case of the low deacceleration parameter qo = 
flo/2 — Xq {Qq, Xq are the density parameter and the cosmological constant, respectively) an average 
flux from sources in a clumpy universe is equal to the flux in the Friedmann-Lemaitre universe 
(flux conservation). For a more general value of 0,q, some authors (e.g., Ehlers & Schneider 1986; 
Peacock 1986) discussed the gravitational magnification probability function by assuming the flux 
conservation. 

Recent observations on high-redshift Type la supernovae (Perlmutter et al. 1999) and on the 
cosmic microwave background (CMB, Spergel et al. 2003) suggest that our universe is accelerating 
in expansion rate. In such a situation, the deacceleration parameter may be no longer small. Then 
it is needed to consider the light ray propagation in a more general inhomogeneous universe model 
with arbitrary values of 0,q and Xq for sources with high-redshifts. 
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In a case with an arbitrary Qq-, we have to take magnification effects by multiple lenses into 
account. In studies of this problem the multi-plane lens theory has been used both in analytical 
approximations (Peacock 1986; Isaacson &: Canizares 1989; Schneider & Weiss 1988a; Wu 1990; 
Marchandon &; Nottale 1991; Seitz &; Schneider 1994) and in numerical simulations (Refsdal 1970; 
Schneider & Weiss 1988b; Watanabe &; Tomita 1990; Rauch 1990; Lee, Babul, Kofman &; Kaiser 
1997; Premadi, Martel, Matzner &; Futamase 2001). In analytical studies many authors (Vietri & 
Ostriker 1983; Pci 1993; Schneider 1993) have assumed that the total magnification by lenses can 
be approximately given by a product of the magnifications of individual lens, and have obtained 
statistically the total magnification by gravitational lenses distributed at random in the universe. 

In this paper we consider the total gravitational magnification factor averaged with all con- 
figurations of lenses distributed at random in a locally inhomogeneous universe and discuss the 
continuous limit (the number A^l of lens planes oo) in which lens planes approach to be continu- 
ously distributed. In §2 the multi-plane lens theory is briefly reviewed and an average magnification 
matrix is obtained in §3. In §4 the continuous limit of the magnification matrix is considered. We 
show that in this limit the average magnification factor satisfies a second order differential equation 

— 1 /2 

and that a angular diameter distance multiplied by {^i{z)) ' {{fj,{z)): average magnification factor) 
reduces to the angular diameter distance of the homogeneous universe (the Friedmann-Lemaitre 
universe) . 



We will give a brief review of the multi-plane lens equation in this section. Suppose that A'l 
lenses are randomly distributed at redshifts Zi{0 < zi < Z2 - ■ ■ < zn^), and that zs = -^a^l+i > -^iVL 
is a redshift of the source (see Fig. 1). The multi-plane lens equation for the source is given by 



where yg denotes the position vector of the source at the source plane and yi is the position vector 
of the light ray at the i-th lens plane (Schneider et al. 1992). In equation (1) D{zi; Zj) is the angular 
diameter distance from the i-th lens plane at Zi to the j-th lens plane at zj. The deflection angle 
OLiiUi) a-t the i-th lens plane is given by 



where Sj(y9 is a surface mass density of the i-th lens and Si denotes the observed region on the 
z-th lens plane. 

We should notice that an image at y^ on the i-th lens plane could be regarded as a "source" 
by the foreground lenses. Therefore the multi-plane lens equation for the "source" at y^ can be 



2. Multi-plane lens equation 




(1) 




(2) 
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Observer Plane 1st lens plane 2ndlens plane A^-th lens plane Source plane 




Fig. 1. — Geometry of multi-plane lens system: Lj is the i-th. lens. Each lens plane is perpendicular 
to observer's line of sight. The origin of each lens plane is set on the line of sight. The light ray 
observed at Uq in the observer plane crosses at in the i-th lens plane. Dij denotes the angular 
diameter distance from the i-ih. lens plane to the j-ih. lens plane (the observer is in the the 0-th 
lens plane and the source is in the {N + l)-th lens plane). 



rewritten as follows: 

= W^f' - (3) 

In the following wc use new variables 9i = yjD{0; zi) (i = 1, • • • , A'^l + 1 = S) which denote the 
angular position of the light ray in the i-th lens plane. Using a dimensionless angular diameter 
distance d{zi\Zj) = D(zi; Zj)/{c/Ho), we can rewrite equation (3) as follows: 

^^ = ^i-E^^«^-[^(0;-^)^^-]- (4) 

Using the x-function introduced by Schneider et al. (see equation [A8] in Appendix A), 
the distance d{zj;Zi) from the j-ih. lens plane to the i-th lens plane is given by d{zj;zi) = (1 -|- 
Zj)d{0; Zj)d(0; Zi){xj — Xi), then equation (4) is rewritten in the following expression 

•i-i 

Gi = 6>i_i - ixi-i - Xi) 1^(1 + Zj)ccj{ej), (5) 

where 

aji9j) = d{0;zj)ocj[D{0;zj)9j] = ^^d\0; zj) J J^d^O'^i: j[D {0; zj)e'j]^ In {Oj - O'j], (6) 
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and O'j and V denote the angular coordinate on the j-ih. lens plane and the observed region, re- 
spectively. An expression similar to equation (5) has been given by Pettcrs, Lcvinc & Wambsganss 
(2001). This recurrence formula (5) determines iteratively the "source" position di in terms of 
9j{j < i) and is useful in a numerical experiment based on the ray tracing method. An equivalent 
equation to equation (1) can be also obtained from the Fermat principle (Blandford &; Narayan 
1986; Kovner 1987): 

Xi,i+i {0i+i - 9i) = Xi-i,i {di - Oi-i) - (1 + Zi)dLi{ei), (7) 

where Xi,j = [Xi-Xj\~^- 

Now we give the magnification matrix ^SjATl in the case of the multi-plane lensing by using 
equation (1) and (4) as 

Ct'l ~[ Xi,S 

where / is the 2x2 unit matrix and Ui and Ai are matrices defined by 

We should notice that Ui defined in terms of cti is slightly different from the matrix Ui (= dcxi/dOi) 
in Schneider et al. (1992). While Ui in their definition depends on the redshift of the source, Ui 
in our definition does not since dj is independent of za,. By virtue of equations (4), (8) and (9), we 
obtain the following form: 

fe=l n=li2=l ifc=l XikM-i--- Ki2mXii,io 
where iq (= ATl -|- 1) > zi > Z2 > • • • > ik-i > «fe > 1, and the matrix Ui can be expressed as 

^'^^'^ ^ ^ci2(O;z,)j^d20'S,[D(O;zO0']t//(0i-0')> (H) 

and 

ri(^) = ^^' r,(^) = ^- (13) 

Recurrence formulae of the magnification matrix given by equation (5) or (7) are also written 



as: 



Ai = - ^(1 + zj)UjAj, (14) 

Xi,i+i {Ai+i - Ai) = xi-i,i {Ai - A-i) - (1 + Zi)UiAi. (15) 
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3. Average Magnification Matrix 



In this section the universe is assumed to be a locahy inhomogeneous, on-average homogeneous 
and isotropic universe in which a mass fraction a (smoothness parameter) of the mean matter 
density p{z) is smoothly distributed, while a fraction (1 — a)p{z) is concentrated into clumps 
distributed at random. The angular diameter distance D{z; z') of this universe from a redshift z to 
another redshift z' satisfies the Dyer-Roeder equation (A6) with < a < 1 (Dyer &: Roeder 1973). 

In this universe a light ray passes through the space with the smoothly distributed mass density 
ap{z) and is gravitationally affected several times by clumps (lenses) located near the light path, in 
general. Since the gravitational magnification factor for the light ray depends on the distribution 
of lenses near the light path, we cannot discuss the individual gravitational magnification factor 
of a source without the knowledge about the configuration of lenses near the light ray from the 
source. Nevertheless it is meaningful to estimate an average gravitational magnification factor for 
light rays which travel in various regions of the inhomogeneous universe, because the factor plays 
an important role in the theoretical analysis of observed data such as m — z relation. 

In the following we consider only the gravitational magnification factor {p) = det {As^n^)~^ 
averaged over all distributions {^i, ■ ■ ■ , ^jvlI of lenses on each lens plane {^f center of the i-th. lens) 
in the locally inhomogeneous universe defined by 



where Q is the solid angle of the observed region T>. Here it should be noticed that to take all 
configurations of lenses into account means to consider observed regions in various directions as well 
as various lens distributions on the individual lens plane. By virtue of equation (10) the average 
magnification matrix is given by 



where Ui{^^]Oi) denotes the matrix for the i-th lens centered on given by equation (11). 

As shown in Appendix B, if the lenses are distributed at random in the universe, i.e., if they 
does not correlate each other, the average of the product of matrices Ui reduces to the product of 
the average matrices \Ui), i.e.. 



(18) 




(16) 




(17) 
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We shall put 9'^ = ^^+0i in equation (11) and assume that all lenses have the same mass profile, 
then all lenses have the same surface density T,{Oi] z-i) = T,[D{0; Zi)0^] which does not depend on 
the center of the i-th lens. Under this assumption, equation (11) becomes to 

AG 



(19) 



and then the average matrix (iJi) can be written by 



Ui{ii,ei)) = 



cHqQ jjj) JJT> 



(20) 



When the matrix Ul{6i — — Oi) is integrated with in a large region, we have \ U[{^j^\6i) 



ah 



T^^ahlQ {a-,h=\ or 2) since the shear terms Fi and r2 in JJ[ vanish because of symmetry. Then we 
find 

-- -^d\0;z^) JJ^d%t{e^;z^)I. (21) 



In equation (21) T,{6i;zi) can be expressed in terms of the matter density p[D{{)\ Zi)6^, Zi] 
in the universe, where is a coordinate along the line of sight given by the cosmological time 
T{zi) and its present value r(0) as c[r(0) — T{zi)\. Since the smoothly distributed matter does 

not contribute to the deflection angle a in equation (6), we can find that the contribution to the 
magnification matrix comes from the inhomogeneous part of p\D{jd\ Zi)0'^, Zj\. Then the surface 
mass density of the i-ih lens plane are expressed as 



T.[D{Q-Zi)e'i] = Sap[D{0;zi)e'i,Zi] ■ \cATi\ = —SaP[D{0; Zi)e'„ Z^] ' 

-Ho [i + Zi)Y [Zi 



where 



6ap[D{0; Zi)el Zi] = p[D{Q- z,)e[, Z^] - ap{z, 
Then we find the mass on the i-th lens plane is given by 



V j^o[i- + Zi)y [Zi) 



Using p{z) = (1 + z)^po and po = SH^rto/STrG, equation (21) can be rewritten as 

U^{^^■,e,)) = -no{l-a) . 



(22) 
(23) 

(24) 
(25) 



Substituting equation (25) into equation (18) and using equations (A8) and (A9), we have the 
average magnification matrix as follows: 



k=l ii=l 12=1 *fc=l 



(26) 
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where 

Ar,, = -Oo(l - «) d{0;z,)Y{z,) ' ^'^^ 

which is the optical depth from the i-th lens plane to the j-th lens plane (i < j). 



4. Continuous limit 

Keeping the total mass of lenses in the universe up to the redshift zs to be constant, we consider 
the limit of TVl ^ oo. In the case of the infinite number of lenses the redshift interval Azi from 
the i-th lens to the {i + l)-th lens plane becomes to be infinitesimal and then the lens plane are 
distributed continuously up to the redshift zs ■ 

Thus, in this continuous limit, summations with respect to is in equation (26) become to 
integrations with respect to ZiS, respectively, and the average magnification matrix is found to be 
given by 

(As)= lim {As,N,)=B{zs)I, (28) 



where the function B{zs) is defined as 



Bizs) = l + J]|--J^o(l-a)| 



' (i + Ci)'d(0;Ci)^i(Ci;Co) 



^ /-^S. (i + C2)^c;(0;C2)rf(C2;Ci) f^"-' {i + CkMo;Ck)d{Ck;Ck-i) 

Jo d{0;a)YiC2) "Jo diO;a.,)Y{Ck) ' ^ ' 

and ^0 = • Equation (29) can be rewritten in form of the integral equation 

o/ ^ 1 3 p (l + ^)2j(0;z)d(^;^g) 

Bizs) = l--no{l-a)j^ d(0^^^)7(^^ 

Prom equation (30) it can be shown that B{z) satisfy the differential equation 

^ |(1 + zmO; z)Y{z)^^B{z)^ + ^Qoil - «) ^ f^^^^"' ^(^) = 0, (31) 
with the initial conditions 

Biz) =1, ^B{z)\ =0. (32) 
z=o az I z=o 

The differential equation (31) can also be obtained by taking the continuous limit of equation (15). 

Since the average magnification factor {^{z)) is given by {niz)) = B~'^{z), now we define a 
new angular diameter distance d{0; z) from the observer to a source at z in terms of (i(0; z) and 
B{z) as 

d{Q; z) = (/x(z))"^/^ d(0; z) = B{z)d{0; z), (33) 
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which is the angular diameter distance magnified with the gravitational lens effect. Prom equations 
(31), (32) and (A6), it follows that d{0;z) satisfies the following differential equation 



d 



dz 



and boundary conditions 



d 



-{{l + zyY{z)-d{0-z)}+-no 



dz 



Y{z) 



d{z) = 0, 



d(0; z) 



z=0 



= 0, 



1. 



z=0 



(34) 



(35) 



Equations (34) and (35) are the same as equation (A6) with a = I. Thus we showed that 
the newly defined angular-diameter distance d(0; z) is equivalent to the angular diameter distance 
dFi^{0;z) in the Friedmann-Lemaitre universe with the density parameter r^o in which all matter 
density is smoothly distributed. 



5. Discussion and Conclusion 

We have to notice that equations (29) - (31), (34) and (35) hold for arbitrary values of fio and 
of a. In the case of the universe with CIq ^ 1, however, the right hand side of equation (29) can 
be understood as the expansion into power series of JIq- The second term Bi{z) of the expansion 
is given by 

R(z)- ^nn r (i + Ci)Mo;Ci)^(Ci;^) ,^ 

Bi{z) - --J^o(l - a) d{0;z)Y{Ci) ^^^^ 

which gives the gravitational magnification effect caused by one deflection. It is interesting that 
—Bi{z) is identical to the optical depth t{z) introduced by Vietri Sz Ostriker (1983). In the case 
of Ho <C 1 it is sufficient to take ^1(2;) into account in order to obtain d{0; z), which is the result 
discussed by Weinberg (with a = 0). 

In a general case of the universe with arbitrary J7o and a, we have to consider B{z) itself which 
includes gravitational magnification effects caused by multiple deflections. The third term -62(2;) 
in the right hand side of equation (29), for example, is written by 



B2iz) 



(1 + Ci)^/(Ci; 4^^^ 1"^' (1 + (,2)^^/(0: C2)rf(C2;Ci) ^^^^ ^3^^ 







d{0;z)Y{Ci) ^'Jo d{0;Ci)Y{C2) 



which is not equal to [Bi{z)]^ /2. In the same manner the (n + l)-th term Bn{z) is found not to be 
equal [Bi{z)\^ jnl. This comes from the fact the total magnification by the multiple deflections can 
not be given by a product of the magnifications by individual deflectors. We have to notice that 
our average magnification factor coincides neither with [1 — t{z)\~'^ given by Young (1981) nor 
with e^'^^^-' obtained by Pei (1993). These differences, however, are not significant in the range with 
2; < 1, but become not to be negligible in the range with z > 1 even in the case of r^o < 1 (see Fig. 
2). 
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Since our universe is locally inhomogeneous, on-average homogeneous, it is needed to know how 
a light ray propagate in the universe. Unfortunately we have no such cosmological model derived 
from the Einstein equation, then we have to investigate which working model to be plausible 
is reasonable and useful to discuss the problem. In this point of view, our conclusion that the 
gravitationally magnified angular diameter distance d{0;z) reduces to dFi,{0;z) in the continuous 
limit is the important result which guarantees the fact that the hypothetical clumpy universe taken 
the gravitational lens effects into account is the reasonable working model to study the light ray 
propagation in the inhomogeneous universe*^. 

This research was partially supported by the Grant-in- Aid for Scientific Research on Priority 
Areas (13135218) of the Ministry of Education, Science, Sports and Culture of Japan. 



A. Background Universe Model and Angular Diameter Distance 

The metric of the Friedmann-Lemaitrc universe (the FL universe) is given by 

els'" = c^dT^ - (T) ( + + sin^ 9d4>''^ , (Al) 

where T is the cosmological time, and a{T) is the expansion factor which has the dimension of 
distance. The Einstein equation in this geometry yields the following relation 

dV SttG , , k^ kc'^ 

where p{z) and A are the mean density of the universe at redshift z and the cosmological constant, 

respectively, and the dot denotes the derivative with T. Since p{z) and a{T) are given in terms of 
their present values pQ and ag by p{z) = po{l + z)^ and a{T) = ao/(l + z), respectively, equation 
(A2) is rewritten as 

Q ' = Hi [Qo(l + zf + Ao - K{1 + zf] ^ HlY\z), (A3) 

where Q,q = SttGpq/SHq, Aq = A(^/3Hq,K = kc^ / HQal{= f^o + Aq — 1) and Hq is the present 
Hubble constant. 

It follows from equation (A3) that the relation between T and z is expressed as 

-cdT = —- (A4) 

Ho (1 + z)Y{z) ^ ' 



''Schneider et al. (1992) have discussed in their book the continuous limit of the multi-plane lens equation with 
the negative surface mass densities and showed the Dyer-Roeder angular diameter distance can be derived from their 
model. 
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Furthermore the relation between z and an affine parameter v along a light ray is derived from 
equation (A4) and from the geodesic equation for the light ray as follows: 

The dimensionless angular diameter distance d{za\ Zh) from a lens at rcdshift Za to another at Zh of 
the clumpy universe in which a mass fraction a of the mean matter density is smoothly distributed 
satisfies the following equation 

(1 + z^fY{zb)j- I (1 + ZkfY{zi,)j-d{za; + ^noa(l + Zkfd{za; z^) = 0, (A6) 

and the following initial condition 

diza;zb) =0, -^d{za;zb) = tt, — \v, (A7) 

z,,=Za [L + Za)Y{Za) 



dzb 



(Dyer & Roeder 1973). The second condition is the Hubble law at redshift Za (Schneider et al. 
1992). 

Schneider et al. (1992) define the x-function in order to express a time delay function by 

{l+Za)d{0;Za)d{0;Zb) 1 

Xab = -77 ^ ^ ■ 

d{Za;Zb) Xa-Xb 

The relation between the x-function and redshift z is also rewritten as 
(see also Seitz, Schneider Sz Ehlers 1994). 



B. Proof of equation (18) 



In this appendix, we prove that the equation (18) holds. In equation (19), the matrix C/j is 
a function of Oi — — Oi. Then, by transforming variables to 0j = + 0j — Oi and putting 
^i(0J = Ui{^i;ei), we have 



II/*"- 11/" 



(Bl) 



where is 5(^j^ ■ ■ ■ ^ij^.)/5(0ii • • • (t>i^) is the Jacobian matrix. We define a 2 x 2 sub-Jacobian matrix 



Jl as 



J ju 
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then the inverse Jacobian matrix can be written as follows: 



7*1 



21 



Til 



4/ 



(B2) 



From equation (4) it can be shown that 0i does not depend on lor k > i and then that 
jh _ p2 

- ■''i2 ~ "' 

form given by 



J*^ = J*^ = • • • = J**' = I and Jj' = O for > ih- Thus we find the inverse Jacobian matrix has a 



O 

\ o 



Jl2 



J^k 

'^ik-1 



o 



(B3) 



I I 



This matrix is a upper triangle matrix of which diagonal component is 1 and then the determinant 
of the matrix is unity. And therefore the determinant of the inverse matrix which is the Jacobian 
of mapping — > 0^ = — — f]^ is unity, too. Thus the right hand side of equation (Bl) can be 
rewritten as 



5(0ii • • • 



From (Bl) and (B4) we have equation (18). 



Vn{(l>i,)---Vi,i4>ik) 



(B4) 
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Fig. 2. — Average magnification factor (/i): Figure (a) shows the average magnification factor in 
the model of = 1.0, Aq = 0.0. Figure (b) does in the model of = 0.3, Aq = 0.7. The thick solid 
line, thick long-dashed line, dashed line and dotted line line show B~'^{z), (1 + Bi)^^ = (1 ~ 
(1 + Bi + S2)~2 and e-2Si(2) = g^"^, respectively. 



